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ABSTRACT 


A new approach to numerical analysis is presented. In- 
stead of providing conventional data one essentially pre- 
scribes the desired pattern and the analysis is performed 
through proportional differences. 

If satellite or other data indicate a synoptic system 
in a numerical analysis should be reanalyzed, a three- 
dimensional manual analysis is performed. Particular at- 
tention is given to the adjustment of the location and 
pee a Lies of the synoptic pattern. No changes are per- 
mitted on the borders of the finite area containing the 
synoptic system. 

Orientation axes for the patterns of the numerical and 
manual analyses are defined by sets of key points. Using 
the space difference between these orientation axes, the 
boundaries and the amount of development specified at the 
key points this scheme guarantees a reanalysis which 
smoothly blends with the original analysis and retains the 


new location and orientation of the synoptic system. 
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I. INTRODUCTION 


Satellites have been providing ever increasingly useful 
data to the meteorologist for synoptic scale analysis of 
weather patterns. Some of the numerous studies which intro- 
duce the terminology and interpretations of synoptic-scale 
cloud patterns as viewed from satellites are references [1], 
fei, {[3), [4], (5], and [6] for video photographs and rer- 
erences [1] and [7] for infrared imagery interpretation. 
However, the major problem has been to determine how satel- 
lite observations can be assimilated with conventional data 
in a numerical analysis scheme, be it infrared or video in 
photographic format or Satellite Infrared Spectrometer 
(SIRS) data. SIRS aeiea is not a big problem in numerical 
analysis programs because the derived heights and tempera- 
tures are essentially conventional data, as indicated in 
[8]. The problems of numerical application of video and 
infrared data from satellite observations are explained in 
reference [9]. One of the earliest attempts to numerically 
integrate the satellite data with conventional data was to 
do it indirectly, through the use of the relative vorticity 
field, as given in [10]. A more objective way of relating 
the synoptic-scale cloud pattern to the relative vorticity 
field is described in [11], [12], and [13]. This relative 
vorticity approach is fundamentally that of Fj¢grtoft [14] 
where the total height field is the sum of the mean-height 


field and the disturbance or relative vorticity field. 





Another early approach was to estimate wind direction 
from the "blow off" of the cirrus plumes of cumulonimbus 
clouds as observed in the video data of TIROS and ESSA 
satellites. The spin-scan cloud camera of the Applications 
Technology Satellites (ATS) can provide video information 
as often as every 20 minutes. The ATS video data may be 
converted into film loops to provide a means of estimating 
wind speed and direction from cloud motion, as explained in 
PES jeand: [Lé}. 

Presently the operational use of inferred data from 
video and infrared imagery requires the use of "bogus" data 
in the analysis program. Bogus data is the name often 
applied to satellite and other non-conventionally observed 
data when these values have been coded in the format of 
conventional synoptic reports for use in a numerical analy- 
Sis program. Most of these bogus data are taken from points 
within an area which has been manually reanalyzed in order 
to provide a more reasonable relationship between the mass 
field and the observed cloud patterns. The bogus data are 
combined with the conventional and SIRS data for a numerical 
reanalysis of the entire area, usually the whole Northern 
Hemisphere, by the standard operational analysis program. 

Mantei and Workman [17] found that preparation of numer- 
ous bogus data values is a laborious process and that the 
determination of the location of bogus data eeuiiae (when 
working with a fixed maximum number of bogus data soundings) 


usually requires several relocations and recomputations to 


? 





finally arrive at the desired synoptic scaled patterns 
of isolines when using a conventional numerical analysis 
program. 

Another solution to the ee of numerical reanalysis 
in sparse data areas iS a numerical scheme based on propor- 
tional differences suggested by Hamilton [18]. In this 
scheme a finite area is delineated for reanalysis with no 
changes in the initial values of the parameters permitted 
on the borders. Within this finite area a meteorologist 
has performed a manual reanalysis at the levels required 
(such as, 1000 mb, 500 mb and 300/200 mb) to establish the 
best spacial and time continuity of all the data available 
within the area. The information available is in the form 
of conventional data (ship, rawinsonde, aircraft, etc.), 
inferred satellite data (synoptic interpretations of satel- 
lite data) and derived satellite data (solutions of sta- 
tistical regression equations for heights of a constant 
pressure surface [13], and heights and temperatures from 
SIRS data, and wind velocities from ATS data). 

The basis of this study is the fact that the most reli- 
able information which can be inferred from satellite imag- 
ery is the location of synoptic scaled systems, such as 
@yclofies, troughs, ridges, jet streams and anticyclones. 
Based on the location of these synoptic scaled systems and 
derived satellite data, manual improvements of numerical 


analyses may be easily accomplished in sparse data areas. 





As indicated earlier, this may be"eequired yy ire wtemeema. 
ent levels within a given finite area of the numerical grid. 

These manual improvements may be incorporated numerically 
when reanalysis is accomplished oe proportional differ- 
ences. The only information required by the computer is the 
location of the finite area, the location of key points with- 
in the initial numerical analysis and within the manual re- 
analysis, and the amount of parameter (height, temperature, 
etc.) changes desired at these key points. Key points are 
those minimum number of points which delineate the synoptic 
pattern within these two analyses. For example, if the 
Synoptic feature to be relocated is a trough, the key points 
locate the axis of the trough (as simulated by straight line 
segments) within the same finite area of both the initial 
numerical analysis and the manual reanalysis. If the synop- 
tic feature is a cyclone, the key points locate the center 
and possibly delineate two axes for pattern orientation. 

The method of this analysis scheme is to relocate the 
values of grid points of the initial numerical analysis based 
on the difference in space between the axes described by key 
points in the two analyses and proportion this difference 
across the finite area, such that there is no displacement of 
the values on the borders. Then apply a four-point Lagrange 
interpolation scheme to the field of displaced values to es~ 
tablish new values at the grid points of the numerical grid. 


The amount of development specified at one or more of the key 





points is proportioned along the axes and then across the 
entire field. The final grid-point values are the sum of 
the interpolated values and the development values. 

A unique feature of this afeaeeae scheme is that the 
desired location and orientation of synoptic features is 
guaranteed to appear in the numerical reanalysis as deter- 
mined from the inferred satellite data, the derived satel- 


lite data and all the available conventional data. 
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IIT. PROPORTIONAL DIFFERENCES SCHEME 


A. BACKGROUND 

As suggested in the introduction, a synoptic-scaled 
System within a sparse data area may require reanalysis 
based on conventional data, inferred data from satellite 
observations, satellite-derived data or any combination of 
these three for one or more levels. A manual reanalysis 
is performed at the required levels to establish the best 
agreement with all the available information, and to ensure 
vertical and time continuity. 

The boundaries of the reanalyzed synoptic feature (cy- 
clone, anticyclone, trough, ridge, or jet stream) are de- 
lineated by grid points of the operational numerical grid 
for the numerical reanalysis using proportional differences. 
For ease of computation in this research rectangular bounded 
areas are used, but this is not necessary as long as the 
border is defined with grid points of the punemedl fig lal 

To ensure a smooth blending with the surrounding areas 
the width of the borders are chosen to be one grid length. 
Figure 1 schematically represents a field of interior grid 
point values chosen for reanalysis. The interior values 
are surrounded by a double border of values which are ex- 
tracted from the total field of grid point values of the 
initial numerical analysis. These border grid point values 
are not permitted to change during the numerical reanalysis 


of this finite area. The initial fields may be any synoptic 
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parameter, but in this study the fields used are height and 


temperature on constant pressure surfaces. 


B. ORIENTATION AXES 

Within the bounded grid are one or two axes which are 
defined by points in order to orient and re-orient features 
on the synoptic scale. These axes are in the general di- 
rections of the i and j axes of the numerical grid. If an 
axis of the synoptic feature is more aligned with the i- 
axis than the j-axis of the numerical grid the orientation 
axils will be referred to as an i-directed axis, and similar- 
ly for an axis oriented primarily in the j-axis direction. 
These two axes can probably describe even the most complex 
patterns which retain synoptic-scale features only. Each 
axis is chosen so that it extends from one side of the in- 
terior grid to the opposite side and is described by a 
Minimum number of key points. The only key points which 
are required by this scheme are the two end penne of the 
axis and any points between the end points where the axis 
changes direction. The axes are defined for both the 
initial numerically-analyzed pattern and the manually re- 
analyzed pattern. The axes of the initial numerical analysis 
are referred to as the analyzed axes and the axes of the 
manual reanalysis are referred to as the observed axes in 
this report. 

If a pattern is to be translated in the i-direction of 
the numerical grid, both an analyzed and an Spooner j- 


directed axis must be specified for the translation to be 


Vs 





accomplished numerically. Likewise, an adjustment of the 
analyzed pattern in the j-direction of the numerical grid 
requires both an analyzed and an observed i-oriented axis 
be specified. Translation in an Sones direction of the 
numerical grid requires both an i and j orientation axis 
be defined for the analyzed and the observed patterns. If 
no translation is desired for a portion of the pattern, the 
analyzed and observed axes are coincident. Figure 2 is a 
graphical representation of a simple translation in either 
the 1 or j direction. In the figure the analyzed axis is 
the solid line and the observed axis is the dashed line. 
The arrows indicate the translation is parallel to an axis 


Sethe numerical grid. 


C. PROPORTIONAL DIFFERENCES 

Within the specified bounded area no changes to the 
initial fields are permitted on the borders while adjust- 
ments of the height and temperature fields are SSNS: ac- 
complished in the interior relative to the axes described 
above. This suggests proportioning the magnitude of the 
adjustment along a row (j constant) or a column (i con- 
stant) of grid points of the numerical grid with the maxi- 
mum at the point of greatest adjustment to no change at 
the corresponding border grid points. This proportioning 
process is performed along an i-axis (row) for translation 
in the i-direction and along a j-axis (column) for trans- 


lation in the j-direction. Translation in an oblique 
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direction relative to the numerical grid is carried out 

in two steps, first in the i-direction and then in the j- 
direction, or vice versa, in order to keep the program 
Simple and efficient with respect Ne computer time. Figure 
3 1s a schematic representation of the two step translation 
for motion in an oblique direction. The diagram illustrates 
the intermediate axes which are determined numerically in 
order for the oblique translation to be accomplished in two 
meeps. that is, the program calculates the location o£) the 
key points of the intermediate axes based on the position 

of the key points of both the analyzed and observed axes. 

im this particular scheme the translation in the i-direction 
is always accomplished first in an oblique movement. There- 
fore, the program calculates the location of the key points 
of the intermediate axes such that after an i-axis shift of 
the analyzed key points only a j-axis adjustment is required 
to place the key points of the intermediate axes:-in coinci- 
dence with the key points of the observed axes. This scheme 
treats the key points of the j-directed axis of the interme- 
diate axis as though they are the key points of a j-directed 
observed axis and performs a standard i-axis shift and de- 
termines a new field of values at each grid point of the 
interior portion of the finite area numerical grid. Using 
the key points of the i-directed axis of the intermediate 
axis as though they are the key points of an i-directed 
analyzed axis the program then performs a standard j-axis 


adjustment. A final determination of the grid-point values 
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resulting from an oblique translation is calculated based 
on the new field of values obtained from the previous i- 
axis shift. In Figure 3 the solid lines are the analyzed 
axes, the dotted lines are the intexnadwae axes, and the 
dashed lines are the observed axes. 

‘Prior to numerical reanalysis using proportional dif- 
ferences, a specification is made for each field to be 
modified of the number of axes, the number of key points 
defining these axes, the motion type (i, j or oblique), the 
amount of development desired (to be discussed later), and 
the location of the boundaries and key points. 

For an 1-axis translation, two j-directed axes must be 
calculated by the program from the given key points of the 
initial numerically-analyzed field and the manually re- 
analyzed field. The orientation axes are assumed to be 
composed of straight-line segments between the defining key 
points. Thus the equations of straight-line segments be- 
tween the key points are determined numerically for both 
the analyzed and observed axes. Using these equations, i- 
values are calculated by the computer along both orienta- 
tion axes for each row (integer j-value) of the numerical 
analysis grid. The difference between these calculated 
i-values is the maximum displacement between the analyzed 
and observed axes for each row. The maximum displacement 
for each row is then linearly proportioned along the row 
from the point of maximum displacement to the locations of 


no displacement, that is, the two opposite border points. 
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Figure 4 is a graphical representation of this proportion- 
ing process. In the figure the solid line is the analyzed 
axis and the dashed line the observed axis. Illustrated 

in Figure 4 are the symbolic form of the proportioning equa- 
fron, the amount of translation (d), the meftiniticn ce: een] 
symbols used in the equation and an illustration of these 


symbols. 


me LNTERPOLATION SCHEME FOR TRANSLATED VALUES 

Once the new location of each initial interior grid val- 
ue (height or temperature) is determined a new value for 
each grid point of the numerical grid must be computed. 
When a single shift is performed, this is accomplished with 
the four-point Lagrange interpolation scheme [19] applied 
only to the initial values of the grid points. This ap- 
proach was judged to give the best results after indepen- 
dent testing using only initial values and a mixture of 
initial and interpolated values in linear and three, FOUr, 
and five-point Lagrange interpolation schemes. Figure 5 
is an illustration of the four-point Lagrange interpola- 
tion scheme. The orientation of the grid points presented 
may be either along a row or a column of the numerical 
grid. If an i-axis adjustment is being accomplished, the 
interpolation is stepped along each row with the i-th grid- 
point always between the second and third known values of 
the parameter being interpolated. Once new values are in- 


terpolated for each interior grid point of a row Cr coum 
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the interpolation scheme is applied to the next row or 
column of the grid. In this manner new values of a par- 
ticular parameter, such as height or temperature, are de- 
termined over the entire grid within the double bordering 
grid points. In addition, this approach guarantees a 
smooth blending of the numerically-reanalyzed field with 


the surrounding initial numerically-analyzed field. 


By. INTERPOLATION SCHEMES FOR VALUES OF DEVELOPMENT 
Changes in height or temperature at the grid points 
due to actual development within the sparse conventional 
data area and those changes in height or temperature re- 
quired to adjust the changes occurring at the grid points 
because of the translation of the synoptic system will be 
summed and called developmental changes, as opposed to 
those changes resulting from the translation reguired for 
proper positioning. Therefore, the value of development 
at any given location within the grid is the value of the 
parameter in the observed (manual analysis) field minus 
the new value of the parameter created by translating the 
analyzed (initial numerical analysis) field. Note that 
with this particular numerical reanalysis scheme the value 
of all key points on a Single orientation axis are not 
changed during a simple translation. Thus the value of 
development at each key point on the observed orientation 
axis is the value of the key point in the observed field 


minus the value of the corresponding key point of the 
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analyzed axis in the initial field. With this scheme the 
value of development need be specified only at one or more 
of the key points on the observed orientation axis. 

Developmental changes are determined numerically along 
the observed axes with a linear interpolation scheme applied 
between the values of development given at the key points 
defining the observed axis. No development is permitted on 
the boundaries. If desired, this program will provide in- 
terpolated values at the defining key points when at least 
one value of development is initially specified. 

A field of development values is determined by linearly 
proportioning the values calculated on the orientation axis 
over a row or column to:the borders. That 1S, an i-axis 
shift results in a row-wise linear proportioning from the 
j-directed cs orientation axis to provide a field of 
development values. Likewise, a j-axis shift results ina 
columnwise linear proportioning from the i-directed axis. 

The development field computations are more complex 
when a field is translated in only one direction but has 
two observed axes, one in the i and the other in the j di- 
rection; or a field is translated obliquely, which by de- 
finition requires two observed axes. In such cases the 
values of development along the j-directed observed axis 
awe calculated but are not proportioned into a field duxeng 


the i-axis translation of the program. When proceeding 
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through the j-axis translation of the program a field of 
development values based on the distribution of development 
Changes along the i-directed axis is determined. Figure 6 
schematically illustrates how the development field based 
on the i-directed observed axis is obtained. The dashed 
line represents the j-directed observed axis extending be- 
tween the key points of this axis. Note that for simplic- 
ity the end points of the axes, which are key points, are 
Met illustrated on the borders. 

The previously calculated values of development along 
the j-directed axis will most likely not agree with those 
interpolated from the field of development values based on 
the i-directed axis. However, the development values cal- 
culated on both the j-directed axis and the i-directed axis 
must be retained in the final field of development. There- 
fore, the field of development based on the i-oriented axis 
must be modified to smoothly incorporate the values on the 
j-oriented axis. To accomplish this smooth blending of 
development values specified along both the i and Jj ob- 
served axes, a four-point Lagrange interpolation scheme is 
used to compute the value of development from the field of 
development values based on the i-oriented axis for each 
row value (integer j value) of the j-oriented axis. A dif- 
ference value is computed peEeeen this interpolated value 
and the previously computed value of development along the 
j-directed axis for each row. These difference values are 


then linearly proportioned along each row to the border or 
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to the i-directed observed axis; whichever is closest to 
the j-directed axis. By definition the difference values 
along the i-oriented axis and on the borders must be zero. 
When the process is carried out for all the rows a field 
of development difference values is generated. The devel- 
opment difference field and the development field based on 
the i-directed axis are added together to arrive at the 
final development field. Figure 7 illustrates the propor- 
tioning of the development difference values from the j- 
oriented observed axis to the border or the i-oriented 
observed axis. Note the existence of "shadow zones" be- 
tween the i-oriented axis and the borders. These shadow 
zones are areas of zero difference which is the correct 
value of difference with these areas. 

The final field of numerically reanalyzed values within 
the finite area is obtained by adding the modified field 
based on the translation of the synoptic pattern and the 
field of development changes. This combined approach guar- 
pieees the location and orientatiom of the synoptic feature 
to be as indicated by all the available information, and 
particularly the satellite imagery data. This numerical 
reanalysis program provides a reasonable estimate of the 
values of conventional parameters at each grid point within 


a sparse conventional data area. 
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F. METEOROLOGICAL CONSTRAINTS 

Initial and reanalyzed fields are checked for vertical 
continuity and conformance with dynamic constraints. The 
first part of the continuity check compares the mean vir- 
tual temperature of the layer with the arithmetic mean 
value of the layer using the temperatures at the top and 
bottom of the layer. The mean virtual temperature is de- 
termined with the 500-1000 mb thickness values used in the 
hypsometric equation. 

Observations from 12 November 1971 to 31 March 1972 of 
the daily Oakland, California, soundings have shown that 
the average difference between these two mean temperatures 
meomesout 2.61°C. This difference in the Oakland soundings 
should be comparable with the inequality of these mean tem- 
peratures in a Marine sounding over a mid-latitude oceanic 
Sparse data area. Therefore, a difference of 2.0°K between 
the mean virtual temperature and the arithmetic mean tem- 
perature is considered within tolerance for the first part 
of the vertical continuity check. The second part of the 
vertical continuity check tests for the existence of super- 
adiabatic lapse rates. 

The dynamic constraint check calculates absolute geo- 
strophic vorticity. Any value less than 1.0 x 107° sec”? 
1s considered to indicate Hae the relative vorticity is 
too negative. 

These checks were used to locate errors while testing 


this scheme. The conventional numerical analysis schemes 
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in operational use have their own checks which ensure proper 
vertical consistency and dynamic constraints of the initial 
fields. Prior to using this scheme as presently programmed, 
the key points should be manually checked for hydrostatic 
and dynamic constraints. If the initial numerical analysis 
fields and the key points are dynamically and hydrostatically 
consistent then this proportional scheme should ensure that 
the final reanalyzed fields will remain meteorologically 
constrained. This has been observed to be the case during 
the preliminary testing of this approach to numerical 
reanalysis. 

If for some reason the reanalyzed fields did not pass 
these or other analysis checks, the new grid-point values 
within the bounded area and the initial grid-point values 
outside the bounded area could be recycled through the 
conventional numerical analysis program as observed synop- 


tic values in order to eliminate analysis discrepancies. 
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Lilie RESUEES 


In the initial phase of this research the various fea- 
tures of the program were extenSively evaluated to ensure 
proper performance. Troughs and ridges (ridges are treated 
identically to troughs in this program) were tested first 
to ensure proper movement in both the positive and negative 
i and j directions. The same was done for cyclones and 
anticyclones (again, anticyclones are treated the same way 
as cyclones in this approach to numerical reanalysis). 
After these features could be moved reliably along an i or 
j-directed axis the programming for oblique movement of 
systems was tested. Finally the development subroutine was 
evaluated for proper performance. 

In order to provide a good test of this numerical re- 
analysis scheme the aa eoanie analyzed fields of height 
and temperature, twelve hours apart, over North America and 
North Pacific Ocean are used to establish initial (Numerical 
Analysis) orsnriaa (Manual Reanalysis) fields. There- 
fore, the heights and temperatures within the finite area 
for reanalysis have undergone a twelve-hour change. A man- 
ual reanalysis is then performed to incorporate the major 
twelve-hour changes of the synoptic system and to blend 
these changes into a no-change boundary around the asso- 
ciated finite area. Thus, the changes near the center of 
the reanalysis areas illustrated in Figure 8 through 1/7 


Simulate the amount of change of height or temperature 
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which would probably be encountered in the operational use 
of this program. All isolines illustrated in Figures 8 
through 17 are delineated on the basis of a linear analysis 
computer program uSing a printer for output. The double 
borders of the finite area of this scheme are delineated 

in the figures for ease of orientation and reference. 

The grid-point values of height and temperature are 
interpolated from the analyzed charts of the National 
Weather Servicewe The griidespacing used in this scheme is 
about 381 km at 60N, which is the same Spacing cumeently 
used in the operational programs of the Fleet Numerical 
Weather Central and the National Weather Service. A grid 
scaezecvotelO by 16, about 1,850 by 3,080 n mi, is illustrated 
in Figures 8 through 17. This provides an area of about 
1,440 by 2,670 n mi for reanalysis. The computer time re- 
quired for the reanalysis of the 500 mb and 1000 mb height 
and temperature fields shown in Figures 10 through 17 is 
about 1 minute on an IBM 360-67, with a total core storage 
of about 175K. 

To further the understanding of this method of numerical 
reanalysis, Figure 8 will be discussed in a step-by-step 
manner. The continuous contours, in this case the isohypses 
of the 300 mb map, represent the initial numerical analysis 
Over a given area. Based eh aeraueenad data, probably from 
a satellite, a meteorologist working in a quality control 
section determines that this initial analysis is in error 


and requires reanalysis. As a first step, he performs a 
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Manual reanalysis at the levels required for spacial con- 
tinuity and establishes the boundaries of the finite area 
of numerical reanalysis. The dashed contours of Figure 8 
simulate this manual reanalysis, but in this specific ex- 
ample these represent the twelve-hour change of the 300 mb 
surface, near the center of the figure, as determined by 
the succeeding National Meteorological Center's operational 
numerical analysis over North America. 

The second step in using this numerical reanalysis 
scheme is the determination of orientation axes, defined by 
the choice of key points, which may or may not be trough 
and ridge axes. In Figure 8 neither of the orientation axes 
are trough or ridge axes. 

It is noted that these orientation axes are i-directed 
and/or j-directed axes and must start on one border and 
terminate on the opposite border. The end key points on 
the borders must be the same for the initial numerical anal- 
ysis and the manual reanalysis (observed field) if the value 
of the parameter being analyzed is not nearly constant along 
the border between the two end points. Along an orientation 
axis are a chosen number of other key points as determined 
by the complexity of the axis. Thism@emplexity is deter- 
mined by the amount of reorientation required and the degree 
Cf GQradHentecontnol ot However, the number of key 
points on an analysis axis must be the same as that on the 


corresponding observed axis. 
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Figure 8. 300 MB HEIGHTS: NUMERICAL ANALYSIS AND 
MANUAL REANALYSIS 
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It should be noted that in cases of simple translation 
of closed systems only a few key points are required. The 
center and the end points of each axis is generally all 
that are required. In a case of uniform filling or deepen- 
ing of a translated, closed system only the center value of 
development is needed. If only filling or deepening of a 
system is required, and this is generally an evenly distri- 
buted value, then only a central value of development is 
required. 

Any extensive pattern reorientation and modification 
will usually require many key points. These points will 
be reguired as anchors for pattern translations and devel-_ 
opment. A good example of such an extensive modification 
is illustrated in Figures 16 and 17. Fine: control of 
trough locations or gradients may require many key points 
on the axes, depending on the axisS orientation, as con- 
trasted to only a few points for coarse control or rather 
Seeavgqhtetroughs. 

The depiction of the analyzed orientation axes, the 
continuous lines through the lettered key points, and the 
observed orientation axes, the dashed lines through the 
primed lettered key points, of Figure 8, indicate a trans- 
lation in only one direction. In this case the positive 
1 direction and then only between key point B/B' and the 
lower, common border key point. That is, when correspond- 


ing key points on the analyzed and observed orientation 
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axesware coincident’ there can be no®translation. Wider. 
tween two coincident points the orientation axes are super- 
imposed upon one another and no translation occurs within 
the data field between these key points. Therefore, the 
remainder of the j-directed axes (B/B', E/E', D/D' and the 
end key point) and all of the i-directed axes (end key 
point, A/A', B/B', C/C' and the other end key point) are 
used for development changes, 1.e., gradient control. Note 
that a key point is not required at the intersection of an 
i-directed and a j-directed orientation axis as might be 
imeessed frem point B/B' in Bigure 6. 

The third step in using this program is the determina- 
tion of the values of development desired. Development is 
specified only at one or more of the observed orientation 
lm=y peints. That is, the primed llettemed points im Figmires 
8, 10, 12, 14, and 16. By definition, development is the 
difference between the value of the parameter after trans- 
lation and the final value in the numerically reanalyzed 
eveld. Since, in theory, simple translation will not clvange 
the value of any key point on the associated analyzed ori- 
entation axis, the value of development of any of these 
key points is simply the difference between the initial 
value and the final value in the manual reanalysis, when 
only one simple translation is performed. Values of devel- 
opment along the other orientation axis, when two observed 
axes are used, or values of development when oblique trans- 


lation is involved must be calculated as the difference 
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between the value determined for the key point after trans- 
lation and the value of the point in the manual reanalysis. 
Since this approach to the problem of numerical reanalysis 
uses a straight forward linear proportion for the changes 
established in the i-axis and j-axis translation, the final 
value of any key point after all translations have been ac- 
complished may be quickly and accurately determined before 
the operation is performed on the computer. Therefore, one 
can always predetermine what this program will produce in 
each step of the computations. This simplicity is not ob- 
tainable when using bogus data in standard analysis programs. 

The values of development used at 300 mb to reposition 
the trough and to properly delineate the axis of the jet 
stream are shown across the bottom of Figure 8. It should 
be remembered that the value of development is always set 
to zero by the program on the inner boundaries of the finite 
area. 

Having determined the location of the finite area and 
the key points along with the desired values of development 
at the selected observed key points, this information 1s 
provided to the program for numerical reanalysis of the 
initial data field. In this case the continuous lines of 
Figure 9 result from the continuous lines of Figure 8. 

The observed orientation axes and the dashed contours 
of the manual reanalysis shown in Figure 8 are reproduced 
in Figure 9 for the convenience of the reader. By com- 


paring the continuous contours with the dashed contours 
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NUMERICAL REANALYSIS ------ MANUAL REANALYSIS 


Figure 9. 300 MB HEIGHTS: NUMERICAL REANALYSIS AND 
MANUAL REANALYSIS 
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of Figure 9 one can evaluate how well the numerical reanal- 
ysis reproduces the essential features of the manual reanal- 
ysis. As can be readily seen, the trough and jet stream 
have been accurately repositioned and delineated without 

the appearance of sub-scale "noise." 

An example of changes at the 500 mb level are illus- 
trated in Figures 10 through 13. The 500 mb height field 
of Figure 10 requires only a slight translation but exten- 
Sive changes in the magnitude of the heights (development). 
As can be seen in Figure 10 the orientation axes are simple 
with only the central key points A and A', and the key 
points E an@ E' required @m the interior for translation. 
The other key points are used solely for development or 
gradient control. The values of development used with the 
cyclone are indicated in the lower right of Figure 10. 
Figure 11 demonstrates the success of this numerical re- 
analysis program with a deepening cutoff cyclone. 

The changes required in the 500 mb temperature field, 
Figure 12, are somewhat more complex than the changes re- 
quired in the 500 mb height field. The center point B, 
key points E, F and the associated primed points on the 
observed orientation axis are requied to adjust the pattern 
through translations while the key points A, C, D, and G 
are required for development and gradient control. The 
actual values of development used in this example are de- 
picted on the right side of Figure 12. The resulting tem- 


perature pattern shown in Figure 11 is close to the desired 
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NUMERICAL ANAENSLS -=->= MANUAL REANALYSIS 


Figure 107 500 Me Berceut.: NUMERICAL ANALYSIS AND 
MANUAL REANALYSIS 


ao 








NUMBRICAL” REANAL Meals =-=-—> MANUAL REANALYSIS 


Figure 11. 500° MB HEhIGhis. ures REANALYSIS AND 
MANUAL REANALYSIS 
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NUMERICAL ANALYSIS ------ MANUAL REANALYSIS 


Figure 12. 500 MB TEMPERATURES: NUMERICAL ANALYSIS 
AND MANUAL REANALYSIS 
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NUMERICAL REAMALYSiIS\---—ae MANUAL REANALYSIS 


Figure 13. 500 MB TEMPERATURES: NUMERICAL REANALYSIS 
AND MANUAL REANALYSIS 
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pattern, but could have been improved with the use of more 
key points. The meteorologist working in quality control 
must decide if the additional information in the sparse 
data area warrants very fine SOnEEOM of the pattern and 
gradients within the pattemayewor aif omly a rough tdeavon 
the overall pattern and gradients are available. He should 
strive not to over-analyze the available data. 

The concluding examples of this numerical reanalysis 
scheme are illustrated in Figures 14 through 17, and are 
on the 1000 mb level. The 1000 mb heights are shown in 
Figure 14. As can be seen in this figure, there is rather 
an extensive pattern change indicated and considerable 
height rises for the 1000 mb pressure surface. It should 
be noted that rather large changes are called for in the 
lower right corner of the grid in view of the proximity of 
the boundary. It is also noted that only three interior 
key points are used to guide the numerical reanalysis pro- 
gram. The values of development are listed in the lower 
left of Figure 14. Considering the many facets of this 
example the reanalysis scheme does surprisingly well, as 
indicated in Figure 15. 

The modification of the 1000 mb temperature pattern, 
Figure 16, is the most difficult numerical reanalysis prob- 
lem of the series, requiring extensive changes to the over- 
all pattern and the creation of a closed center of cold air. 
The choice of key points in this case clearly illustrates 


the fact that key points need not be on an actual trough or 
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NUMERICAL ANALYSIS ------ MANUAL REANALYSIS 


Figure 14. 1000 MB HEIGHTS: NUMERICAL ANALYSIS AND 
MANUAL REANALYSIS 


44 








NUMERICAL REANALYSIS ------ MANUAL REANALYSIS 


Figure 15. 1000 MB HEIGHTS: NUMERICAL REANALYSIS AND 
MANUAL REANALYSIS 
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ridge axis, and in some cases cannot be if the desired mod- 
ifications are to be achieved. 

The extreme protrusion of warm air required in the mod- 
ification of the upper peril coenmor tne figure could not have 
been accomplished with development changes alone. Note 
that the i-axis. translation is both positive and negative 
(G to G'). The required development is given on the right 
Side of Figure 16. The results in Figure 17 show how well 
the above modifications worked. The extreme protrusion of 
Warieait cequired in the Geen geet Lepneand the mmauked ceolind 
between points C and D resulted in the corollary develop- 
ment of a small warm center because of the proportional 
manner in which this program spreads changes in the data 
field. However, the difference between the sharp protrusion 
and the closed center would probably not be significant in 
an operational analysis. 

The effect of a mistake in entering values of develop- 
ment can be seen in Figure 17. The value of development 
mequared for key point I' as’ -3°. The@manus sign was in- 
advertently left off so the value of development actually 
MSedqd in the program is +3°. Gas resultegameehe slight 
Eldging effect" in the vicintiy wer 1' ee igure 17. 

This series of figures illustrates that the degree of 
accuracy desired is determined by the number of key points 
chosen to establish the changes to be effected ee ie 


program. The complexity of the reorientation of the pattern 
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NUMERICAL ANALYSia: ---=— MANUAL REANALYSIS 


Figure 16. 1000 MB TEMPERATURES: NUMERICAL ANALYSIS 
AND MANUAL REANALYSIS 


47 








NUMERICAL REANALYSIS ————.. MANUAL REANALYSIS 


Figure 17. 1000 MB TEMPERATURES: NUMERICAL REANALYSIS 
AND MANUAL REANALYSIS 
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and orientation of the various axes also bears heavily on 
the number of key points used. These examples have, how- 
ever, demonstrated that the program is flexible and yet 


simple and easy to use. 
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iV. CONGEVUS fells 


A significant aspect of this methodsas that the spocaemen 
and orientation of cyclones, anticyclones, troughs, ridges 
and jet streams is guaranteed in this method of numerical 
reanalysis, and no such statement can be made with the oper- 
ational use of bogus data in the standard analysis programs. 
The information which is directly interpretable from video 
and infrared imagery is the position and orientation of cy- 
clones, troughs, ridges, jet streams and anticyclones, and 
not quantitative values of bogus data at the numerous points 
required in conventional analysis programs. This reanalysis 
scheme concentrates on the information which can be most 
reliably taken from this type of satellite data and minimizes 
the use of bogus data values which must we aneewEwd from 
statistical regression equations, space continuity, time 
Gontinuity and dynamical constraints. 

This scheme appears to be easier and faster to use than 
the current operational procedures using bogus data. In 
mediton, the quality control meteorologist knows exactly 
what the pattern and central or key values will be with this 
method of numerical reanalysis. 

Based upon this preliminary research, this numerical re- 
analysis scheme appears to offer an operational method of 
improving the "update analysis", the "first guess field” or 
both of these fields, until quantitative soundings of the 


atmosphere are routinely available over a dense network of 
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data points from satellite-borne anstruments. The "update 
analysis" is the last analysis performed by a numerical 
center for a given synoptic time and, for the upper-air 
analysis, the "first guess field" is essentially a twelve- 


hour prognosis based on the "update analysis‘" 


ak 





Lie 


BIBLIOGRAPHY 


Air Weather Service Technical Report 212, Application 
of Meterological Satellite Data in Analysis and 
Forecasting, by Anderson, R. K., and Others, June 
1969 (supplement, 1971). 


National Environmental Satellite Center, ESSA Direct 


Tiron sitines S 1 ON oy Ss Peluoe hs Gurce, = ooo. 


Naval Air Systems Command Project FAMOS Research Re- 
port (4-67), Guideg@ror Interpretation of Satellite 
Photography and Nephananalysis, August 1967. 


Navy Weather Research Facility Report 33-1169-148, 
The Use of Satellite. Piceumes for Suntace and, 500 
mb Chart Analysis, November 1969. : 


Office of Naval Research, Branch Office London (Eng- 
land) Technical Report ONRL-44-66, Operational 
Guide to Synoptic Applications of Meteorological 
Observations from Satellites, by Hamilton, H. D., 
November 1966. 


World Meteorological Organization Technical Note No. 
75, The Use of Satellite Pictures -in Weather 
Analyses and Forecasting, by Anderson, R. K., 
Herguson, E. Wate cicmOlivon ee Ve J., 1966. 


Navy Weather Research Facility Report F-0970-158, 
Guide for Observing the Environment with Satellite 
infrawed Imagery 2ay Bittner, tease, ama Rugg@hes , 
K. W., september 1970. 


Navy Weather Research Facility Technical Paper No. 
25-67, A Comparison of Maritime Vertical Tempera- 
ture Profiles, Radiosonde Versus Satellite Infrared 
Spectrometer (SERS), by Naisestien, Ro My Cetober 
1969. 


National Environmental Satellite Center, ESSA Tech- 
nical Report NESC-41, The SINAP Problem: Present 
otatus and Futuve Prospeees, bye ecClain, E. P., 
Setrober 1G. 


National Environmental Satellite Genrer Meteorolegreal 
Satellite Laboratory Report 36, Experimental Use of 
Satellite Pictures in Numerical Prediction, 
Broderick, H. J. (’Mec@iaan rhe ond nu zee ee, 

Me A., January tae 


2 





bliss 


a 


3 


14. 


iS. 


ies 


ey 


3 


Le 


Meteorology International, Ime., Final Report, Con— 
- tract No. E-93-67(N), An Approach to the SINAP 
Problem: A Quasi-Objective Method of Incorpora- 
ting Meteorological Satellite Information in 
Numerical Weather Analysis, by Nagle, R. E. and 
Clark, J. Rs) eicsc. 


Project FAMOS Research Report 1-68, Techniques for 
the Objective Assembly of Meteorological Satellite 
Data, by Nagle, peat anc Cilaia a) ee oigio 


National Oceanic and Atmospheric Administration 
(NOAA) Technical Report NESS-55, The Use of 
Satellite-Observed Cloud Patterns in Northern- 
Hemisphere 500 mb Numerical Analysis, by Nagle, 
Ree becand Hayeent ermine, 1O7l. 


Fjg@etort, R.,; “Ona Nomemmcal Method of Imteqrating 
the Barotroplic Vermiemtysnduatten, Eins, Vea 
pe 179-194, JabesAvc Meise ; 


Semebreny, S. M. and Others, Comparison of Cloud 
Motion Vectors and Rawinsonde Data, Stanford 
Research Inst., Menlo Park, Calif., December 
1969. 


Senebmeny,.S. M., Wemement. J2, and-Hadfield, Rote. 
Further Comparisons of Cloud Motion Vectors with 
Rawinsonde Observations, Stanford Research Inst., 
Menlo Park, Calif, 2) Jueme: 1970. 


Mancer, ©. J. and Workman, ©. Es, Beecrimentalyuseme.r 
Satellite Data in Numerical Analysis and the Effect 
On a Primitive Bevatvon Prediceion Scheme, M.S. 
Thesis, Naval Postgraduate School, Monterey, Calif., 
September 1971. 


Office of Naval ReSe@arch weranch Office London (Eng- 
land) Technical Report ONRL-64-77, The Second 


Operational Step Towards Complete Numerical inte- 


gration of Satellite Observations with Conventional 
Data, by Hamilton, H. D., December 1967. 


Todd, John and Others, Survey of Numerical Analysis, 
p. 36, McGraw-Hill, eu 


3) 





INITIAL DISTRIBUTION LIST 


Library ,-¢ede 02 
Naval Postgraduate School 
Monterey, California 93940 


CDRH Pea iani) ton, secteur 
Department of Meteorology 
Naval Postgraduate School 
Monterey, Californma 93940 


Professor Robert J. Renard, Code 51Rd 
Department of Meteorology 

Naval Postgraduate School 

Monterey, California 93940 


LCDR R. B. Glaes 
USS Constellation (CVA-64) 
FPO San Francisco 96601 


Naval Weather Service Command 
Naval Weather Service Headguarters 
Washington Navy Yard 

Washington, D. C. @20390 


Commanding Officer 
Fleet Numerical Weather Central 
Monterey, California 93940 


Officer in Charge 

Envixronmental Prediction Research 
Facility 

Naval Postgraduate School 

Monterey, California 93940 


befense Document at von center 
Cameron Station 
Alexandria, Virginia 22314 


Atmospheric Environmental Service 
Forecast Research Division 

4905 Dufferin Street 

Downsview, Ontario, Canada 

Attn: Dr. F. Bl Muller eect we.) 


54 


No. Copies 


Z 


10 





Oe 


jlo 


Dynamic Prediction Research Unit 
c/o Central Analysis Office 
Department of the Environment 
Box 158, Montreal AMF 
Montreal International Airport 
Dorval, Quebec, Canada 

Attn: Dr. lan) Deeruehemiend 


Commonwealth Meteorology Research Centre 
Mackay House 

380 Russell Street 

Melbourne, Vic., Australia 3001 

Attn: Miss F. J. Taubman 


a 





Secunty Classification 7. 
DOCUMENT CONTROL DATA-R& 


and indexing annotation must be @n 


> 


= , 


tered when the overall report is classilied) 


(Security classification of title, body of abstract 
ZIGINA TING ACTIVITY (Corporates author) 2a. REPORT SECURITY CLASSIFICATION 
aval Postgraduate School Unclascautee 
paeemewe California 93940 25. CRORE 
EPORT TITLE 
1f£ferences 


umerical Reanalysis through Proportional D 


L\ESCRIPTIVE NOTES (Type of report and, inclusive dates) 


laster's Thesis: September 1972 


.U THORS) (Firs( name, midale initial, jest name) 








Roger Burton Glaes 
REPORT OATE 
Dl 


CONTRACT OR GRANT NO. 


PROJECT NO. 


Ul 
° 


1 DISTRIBUTION STATEMENT 


Approved for publie release, 


1. SUPPLEMENTARY NOTES 


Naval P 


Ja. TOTAL NO. OF PAGES 


12. SPONSORING MILITARY ACTIVITY 


7b. NO. OF REFS 


Ike 


September 1972 | 
9a. ORIGINATOR'S REPORT NUMBER(S) 


96. OTHER REPORT NO(S) (Any other numbere that may bo asaigned 
thie report) 





distribution unlimited. 


ostgraduate School 


Monterey, California 93940 


3. ABSTRACT 


A new approach to numerical ana 


viding conventional dat 
and the analysis is per 


lysis is 


Tf satellite or oth 
cal analysis should be reanalyzed, 
is performed. Particular atten Eonmis 
location and orientatil 
mitted on the borders of t 


system. 


A one Essentially prescE 
formed through proportional differences. 


er data indicate a synoptic s 
a three-dimensio 
given to the adjustment of the 


on of the synoptic pattern. 
he fintee area contaLnan 


presented. Tnstead of pro} 
‘bes the desired pattern} 


ystem in a numeri- 
nal manual analysis 
No changes are per~ 

g the synoptic 


he numerical and manual anal- 


Orientation axes for the patterns ONE 
Using the space difference 


points. 


yses are defined by sets Of ey 
the boundaries and the amoun 


between these orientation axes, 


analysis which smoothly 
the new location and orie 


nam Tet 


Ri teraea Pee APNE) . 
e 


S/N 0101-807-6611 


velopment specified at the key points this sche 
blends with the origina 


ntation of the synoptic system. 


t Ofnaes 


me guarantees a ren 
1 analysis and retains 


ecurity Classification a-31408 





Security Classification 


KEY WORDS 





Numerical Weather Analysis 
Meteorological Satellite Observations 
1000 mb, 500 mb, and 300 mb Analysis 


- ~ 2 


DD L473 (Back) 


S/N 0101-807-6821 57) 


LINK ¢ 


= " 


| 
| 
| 


Security Classification 


A-31340>% 


VAST ee ee 


ESOT ANS A RIENCE 











esis 135543 
4565 Glaes 

opm Numerical reanalysis 
through proportional dif- 
ferences, 





thesG4565 


DUDLEY KNOX LIBRARY C. 





